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We found an amazingly simple general ’’all-band” intensity profile of bandgap (Bragg) solitons 
for arbitrary parameters of spatially-periodic nonlinear systems, similar to those of multi-frequency 
stimulated Raman scattering, in particular the so called Lorentzian-proHle solitons. We also found 
nonlinear eigen-modes of such system that propagate without energy exchange between waves. 
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I. INTRODUCTION 

In nonlinear optics, solitons, both temporal and spa¬ 
tial, became familiar objects, which can be originated by 
quite a few fundamental nonlinear processes. The ones 
of interest here are so called gap (or Bragg) solitons [1-5] 
(for review, see e. g. [6]) that are due to the interac¬ 
tion of light with spatially-periodic nonlinear structures, 
whereby counterpropagating waves are strongly coupled 
to each other via distributed back-and-forth Bragg re¬ 
flection. The Bloch theory of linear periodic structures 
predicts the existence of ’’prohibited” zones - bandgaps 
- where due to resonances between incident wavelength 
and the period of spatial modulation, the incident wave 
cannot propagate for long, and the structure becomes 
almost fully reflective; those are essentially multilayered 
Bragg reflectors. The spectral width of such a bandgap is 
proportional to the contrast between refractive indices of 
constituent materials (or waveguide). When these refrac¬ 
tive indexes also depend on the light intensity (e. g. due 
to Kerr-like nonlinearity - self-focusing if positive, self- 
defocusing if negative), the system may exhibit a rich 
host of nonlinear effects, some of the most interesting 
being gap solitons that emerge as a critical phenomenon, 
whereby the light with its frequency being inside the 
bandgap, ’’pushes” the bandgap edge away and carves 
conditions for itself to penetrate deep into the structure 
by forming gap solitons. 

A new property of bandgap solitons, predicted first in 
[1], which gave them a distinct place amongst the soliton 
crowd, was that they may exist as standing (immobile, 
stationary) field objects. In fact, they might be viewed as 
a simplest, ID-self-trapping of light. It was demonstrated 
then [2] that under certain conditions those solitons could 
be similar to familiar nonlinear solitons in regular fibers 
[3] and 2D-self-focusing [4] and be governed by a rescaled 
cubic Schrodinger equation. 

In further development, it has been shown [5] that 
Bragg systems can support even more general gap soli¬ 
tons that have non-zero but very slow group velocity, 
giving rise to nonlinear ’’slow light” objects. The slow 
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bandgap (SBG) solitons have later been observed exper¬ 
imentally in [7], and immobile solitons - in [8]. Gap soli¬ 
tons have by now been shown to emerge in many other 
bandgap systems, such as e. g. in a Bose-Einstein gas [9] 
and near critical point in underdense plasma due to rela¬ 
tivistic nonlinearity of electrons [10]. While SBG-solitons 
have been shown [5] to exist within a linear bandgap, an 
approximate analytical solution in [5] have been found 
only at a close vicinity of the bandgap edge (which in 
the case of self-focusing nonlinearity we call here a ” blue 
edge”, see Section IV below). Its intensity profile was 
again similar to nonlinear Schrodinger equation soliton. 
(As one can see below. Section IV, it is the weakest and 
longest soliton.) A greatly important development in the 
field was a proof [11] (based on generalization of a so 
called ’’massive Thirring model”, see references in [11]) 
that related coupled nonlinear partial differential equa¬ 
tions are fully integrable, which in particular allowed for 
a complex solution for stationary SBG-solitons with an 
arbitrary detuning within a bandgap. 

In this paper we show that the family of intensity pro¬ 
files of all the solitons within the entire Bragg bandgap 
are described by a simple elegant general formula (see 
Eqs. (3.6) and (4.1) below), which may be called an all¬ 
band Bragg soliton, derived by us based on a standard 
Bragg wave equation (Section II) and directly looking for 
a solution that satisfies the conditions of a fixed group 
velocity and bright-soliton asymptotic (Section III). It 
can be shown that these solutions are consistent with 
the general results [11], yet their analytical simplicity 
and transparency related to the use of intensities and 
phases as output variables, allowed for easy visualization 
and analysis of the soliton intensity profiles, in particu¬ 
lar in the less studied, and perhaps the most interesting 
area near red-edge of the bandgap (for a positive Kerr- 
nonlinearity), where they assume an extremely simple yet 
unusual for solitons Lorentzian profile, see (4.5) below, 
which was missing in [11]; they are similar to a limiting 
case of Raman scattering solitons [12]. This approach 
also brings up all the relevant invariants of motion. We 
analyzed in detail both moving (Section IV) and station¬ 
ary solitons (Section V). Using a similar approach for cw, 
non — soliton, solution, we found cw — eigenmodes of the 
system, in particular in a ring Bragg reflector (Section 
VI), and discussed the ramification and possible applica¬ 
tions of our results (Section VII). 
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II. WAVE PROPAGATION MODEL: BRAGG 
REFLECTION + KERR-LIKE NONLINEARITY 

Considering a ID-model (e.g. as in an optical fiber), 
we assume that the periodic grating in it is formed by 
the modulation of the refractive index, n (this can also 
be done by slight periodic tapering of fiber, which can 
be dealt with by a similar math-description.) In a uni¬ 
form, unmodulated (n = const) fiber the electric field E 
propagation is governed by a regular wave equation 

d'^Ejdz'^ - {n^lc^)d‘^Eldf = 0 (1.1) 

which in low-dispersion case can be decomposed into the 
set of two first-order equations for ’’forward”, Ei, and 
’’backward”, E 2 , traveling waves as: 

{-lydEj/dz - {nlc)dEjldt = 0; j = 1,2 (1.2) 

whose w-monochromatic radiation solution is Ej = Aj 
exp[—z(—l)'^fcz — iijjt]/2 + C.C., with Aj = const, and 
k = Lo/cn = 2'Kln\, n is an unperturbed refractive in¬ 
dex, and A is a free space wavelength. (As is common 
in the bandgap theory, we will neglect here the intrinsic 
dispersion of n{uj) in a fiber, since the bandgap disper¬ 
sion due to index modulation greatly exceeds that of a 
regular waveguide [5].) 

A spatially modulated fiber making a periodic grat¬ 
ing with a period Lb, has its resonant Bragg wavenum¬ 
ber, free space wavelength, and frequency respectively 
as ks = t:/Lb, Xb = ^uqLb and ujb = kBc/no] a 
normalized half-width of a linear bandgap around lob 
is fi = An/no <C 1, where An is the index modu¬ 
lation amplitude. No radiation with is al¬ 

lowed then to propagate in a sufficiently long system; 
here A = w/ws — 1 is a normalized Bragg detuning. 
Furthermore, in the presence of a Kerr-like nonlinear¬ 
ity, refractive index depends on the intensity of light as 
n = no -I- n^A^, where n^ is a coefficient due to 
(Kerr) nonlinearity. The total refractive index can be 
written then as 

n = no + Ancos{2kBz) + nxE'^ (1.3) 

Generalizing (1.1) to have the index n as in (1.3), and 
seeking its solution as the sum of traveling waves, E = 
El + E 2 , with 


(nK/no)(|A,f+ 2|A3_,f)A, =0; j = l,2 (1.5) 


where the factor 2 in the nonlinear term is originated by 
intensity-induced non-reciprocity [13]. In linear case Eq. 
(1.5) is similar to coupled equations for counterpropagat- 
ing modes in a Bragg reflector [14]. 

To normalize (1-5), we use the scales for intensity, 
E'/qj^ = An/nK, distance, znl = 1/m^b, and time, 
tNL = I/mi^b- Introducing then dimensionless envelopes, 
Gj = Aj/Enl, time t = t/tNL, and distance C = z/znl, 
we rewrite (1.5) as: 


i 



da-j 

dr 


dj 5 a3_jH- 


SK (Uj + 2 ul_j) Gj =0; = \Gjy (1.6) 

where sk = sign(nK/An). For a fixed sk, only one con¬ 
trolling parameter, a normalized detuning 5 = A/ji = 
(uj/lzb — l)Bo/An, is left in (1.6). A linear band is con¬ 
stituted then by the condition 5'^ <1. 


III. WAVES WITH FIXED GROUP VELOCITY 

As a next step let us find the solution of (1.6) for the 
class of modulated coupled waves with a fixed group ve¬ 
locity, P = Vgrno/c, (/3^ < 1); the solitons will be part 
of that family. To that end, we assume that in the frame 
moving with that velocity, i. e. ’’comoving frame”, the 
fields are time-independent. Introducing then time-space 
variables in that frame as C /3 = ^ ~ Pt, Tg = t, we are 
looking for the solutions with d/drp = 0 to derive ordi¬ 
nary nonlinear differential equations 

— z[/3-t-(— l)yGj-\-SGj-\-Go—j~i~SBl[Uj-\-2u^_j)Gj =0 (2.1) 

where ’’prime” stands for d/dt/fj. To elucidate direct re¬ 
sults for the intensity profiles, from this point, as differ¬ 
ent from [11], we will be using an ” amplitude&phase” 
approach similar to our previous work [15,10], which 
greatly simplifies the calculations. To that end, we write 
Gj = Uj exp(zC'j) with Uj and pj - real, rewrite (2.1) as 

~[P + ~ y/^j) + Uj[5 + sxiu^ + 2.ul_j)] + 


Ej = Aj {t, z) exp[—za;B(l -I- A)t — i{—iy kBz]/2 + c.c.; 

(1-4) 

with j = 1,2, presuming the envelopes Aj to vary slowly 
in time and space, as JAnj, <C no, and neglecting 

higher-harmonics generation (a common approach [1-6] 
in nonlinear Bragg reflection because of large difference 
in phase velocities), we obtain truncated equations of evo¬ 
lution being nonlinear-l-modulated counterparts of (1.2) 
for those envelopes as: 

i [—{—l)^dAj/kBdz + dAj/ujBdt] + AjA + fj.A3-j + 


= 0 ( 2 . 2 ) 

and separate real and imaginary parts in Eq. (2.2): 

[1 -I- {-l)^P]uj' = U 3 -j sin(())i - P 2 ); (2.3) 

— [<5 + sk{u^ + 2u3_j)] = 

-U 3 -j COs(C'l - p 2 )- 


(2.4) 
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Using eq (2.3) for both j = 1 and j = 2, we have {u\ — 
u^y = /3{ui + ^^ 2 )^ hence the invariant of motion for this 
class of solutions: 


{ul — u\) — I5{u\ + u\) = D = inv (2.5) 


Using this in (2.3) and (2.4) and introducing ’’combined” 
variables 


S = uf+U2; ^ = 4>1 — 4’2', U = ^l+^2; P = UiU2 

( 2 . 6 ) 

we see that their spatial dynamics is governed by equa¬ 
tions: 


S'= 


4Psin$ 
1-/32 ’ 




$' = - 


2(,5 + skS) + [^(1 - /32) - /3D]{sk - cos$/P) 


1-/32 


( 2 . 8 ) 


_ 2(3{S + skS)-D{sk-cos^/P) 

S (2.9) 


From eqns (2.7) we obtain phase space equation for P, S 
independent of $ and S, which greatly simplifies further 
calculations by separating those variables: 


dS 4P 

dP ~ 5(1 - /32) - /3P 


( 2 . 10 ) 


whose solution, e. g. for P{S), is immediately found as 


PiS) = V(1 - P^){S - 5™„)[(5 + 5™„)/2 - f3D]/2 

( 2 . 11 ) 

where Smin is an integration constant chosen in such a 
way that P = 0 when S = Smin- Having this result and 
using (2.8) and first eqn (2.7), we obtain now a phase 
space eqns for $,5, as 

d^/dS = (2.12) 


2(,5 + skS) + [5(1 - /3") - PD][sk - cos$/P(5)] 
4P(5) sin $ 


IV. ALL-BAND MOVING SOLITONS 


so that eqn (2.12) for variables 5 and $ reads now as: 

M _ 2(35^5(3-/3^)-2a/1-/32 cos4> 

dS ~ 25sin$A/l-/32 

the solution of which is readily found as cos$ = C/S + 
[(5 -I- 5(3 — /3^)/4]/a/ 1 — /32 where C is yet another inte¬ 
gration constant. Due to (3.1) we set C = 0, so finally 

cos $ = [(5 + sk5(3 - /32)/4]/\/1 - /32 (3.4) 

Using (3.4) and (3.2) in the first eqn (2.7), we obtain a 
first order eqn for one variable, 5(^/3), alone 


5' = 25 


yiHZ 


[(5 + sk5(3-/32)/4]2 
1-/32 


(3.5) 


whose elegantly simple yet little familiar in the soliton 
theory solution can be readily found: 


Spkjl + B) ^ _ Ssk 

cosh[a(C/3 - Co)] + B ’ i/l - /32 


(3.6) 


where combined peak intensity, Spk , and size-related pa¬ 
rameter a are 


Spk — 


4\/r^(i-p) 

3-/32 


2v'l-/32-(52 

“ 1-/32 


(3.7) 


In (3.6) the soliton peak position, Coj is yet another in¬ 
tegration constant; without loss of generality we will as¬ 
sume Co = 0. 

Amazingly, profiles (3.6) that suggest a broad control 
of the shape and width of a soliton determined by pa¬ 
rameters a and B (notice here that B can be either neg¬ 
ative or positive), coincide with the soliton profiles that 
show up in the multi-cascade stimulated Raman scatter¬ 
ing [12], with completely different mechanism of their for¬ 
mation. Apparently, such solitons are common solutions 
for many situation and system with interference between 
waves with different wave-vectors: in those two examples 
either co-propagating waves (as in Raman scattering), or 
counter-propagating waves (as in Bragg scattering). 

Individual waves intensities are now found via (2.5) as 

u? = 5(1 +/3)/2; = 5(1 -/3)/2. (3.8) 


The solution of (2.12) in the form e. g. $(5), being 
substituted in (2.7) and (2.8), after integration yields 5, 
P, and $ as function of (/p. In general, all of them can 
be shown to be periodic functions of the distance (/p, 
except for solitons. The so called bright solitons must 
then satisfy the condition, 

Uj —>■ 0 at Id —>■ 00 (3.1) 

From (2.5) then P = 0, hence uf — + “i)- 

From (2.10) and (2.11) we have also that Smin = 0 and 

(3.2) 


Substituting (3.6) into (3.4), we have an expression for 
the phase difference ^ = <j>i — <j> 2 - 


1 + B cosh(a^/3) 
P-I-cosh(a^^) 


or 


ton 4. = ,3 3 ) 

\/l - /3^[1 -H Pcosh(aC/3)] 

At the peak, t/p = 0, we have $ = 0, so the counter- 
propagating waves are of the same phase at that point, 
(/i = (j) 2 - At (p —>■ 00 , we have tan<l> —>• a/1 — /32 — 6^/S. 


P = 5V1-/3V2, 
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At the exact resonance, <5 = 0, $ is reversed by tt as we 
go from = —oo to = oo. Tedious but mundane 
calculations show that Eqns. (3.7) - (3.9) are consis¬ 
tent with results [11] (see Eq. (6) in [11]); in particular, 
the major soliton parameter Q in [11] is related to the 
parameters used here via cosQ = B. The simplicity 
and transparency of (3.7) - (3.9) however, allows one to 
easily analyze the soliton intensity profiles, in particu¬ 
lar in the less studied, and perhaps the most interesting 
area near red-edge of the bandgap (for a positive Kerr- 
nonlinearity), where they assume an extremely simple yet 
unusual for solitons Lorentzian profile, see (4.5) below. 

Notice from (3.6) - (3.8) that the bandgap for the mov¬ 
ing solitons with a given velocity, /?, is narrower than in 
the linear case: 5^ < 1-/3^. By the same token, near 
both the edges of a bandgap, 1 — <C 1, only slow soli¬ 
tons are allowed, since then 0^ < 1 — S^. On the other 
hand, in the middle of the bandgap, <5 = 0, the solitons 
are allowed with the speed up to the maximum /3^ = 1. 


V. STATIONARY (IMMOBILE) SOLITONS 


For immobile, or stationary, solitons, /3 = 0, we have 
ul = U 2 = Sl2 (hence the analogy to standing waves), 
and the combined intensity profile, S, and phase $ are 
found from (3.6) and (3.9) as 


4(l-(5")/3 
S + cosh(2^-\/l — S^) ’ 


(4.1) 


yr^sinh(2cyr^) 

tan<l> = - . — (4.2) 

1-b (5cosh(2CyF^) 

Since this analytical solution in the entire bandgap, i. e. 
for any < 1 is now available, one can analyze how it 
evolves as the laser frequency is tuned from upper (blue) 
edge of the bandgap, <5 = 1, to the lower (red) edge, 
(5 = — 1. For example, close to the blue edge, 0 < 1 — (5 = 
AS <C 1, we have a low-intensity long pulse (see Fig. 1, 
lower curve): 


4(A,5)/3 


cosh2(^-\/2A(j) ’ 


$ « y^tanh(Cy^) (4.3) 


consistent with [5], which is a familiar soliton of a cubic 
Schrodinger equation, whereas at the exact Bragg reso¬ 
nance, (5 = 0, we have a somewhat different profile. 


C 4/3 
cosh(2^)’ 


tan$ = sinh(2C) 


(4.4) 


It is worth noting that this is a solution of a generalized 
Schrodinger equation with the higher order nonlinearity 
of the 5-th order, whereby n = no + n^E'^ [16]. 
Finally, at the red edge, 0 < 1 -b 5 <C 1, the main body 
of a standing soliton (with the maximum peak intensity) 
has the most unusual, Lorentzian, profile: 


S « (8/3)/(l + 40); tan $ « 4C/(1 + 40) (4.5) 


« 



- 3 - 2-10123 
Normalized distance, ^ 


FIG. 1: Stationary solitons at various detunings 5 (S^ < 1) 


similar to a limiting case of stimulated Raman scatter¬ 
ing solitons [12]. Of course, far away from the peak, 
0 ^ (1 — S^)~^, its profile decays exponentially, S oc 
exp(—2C-\/l — S^). 

The total power carried by a stationary soliton in gen¬ 
eral case, i. e. for arbitrary S, is 

W{6)=j SdC = ^ tan-i (4.6) 

and the soliton width in the C-axis at the half-peak in¬ 
tensity is 


Zl/20) — 


cosh ^(2-b5) 


(4.7) 


At 5 > 0, AS = 1 — (5 <C 1 (the blue edge in case of 
Kerr-nonlinearity), solitons are the weakest and longest: 


5, 


pk 


4(A(5) 


W : 


4VM 


^ ^ cosh ^(1) 

^1/2-(4.8) 


3 ’ 3 ’ V® 

In the middle of linear band, 5 = 0, we have 

Spk = 4/3; IT = 27r/3; Z 1/2 = cosh“^(2) (4.9) 


The solitons at 5 < 0 and I -b 5 <C 1 (red edge) are the 
strongest and shortest ones, albeit not by much: 

Spk^ 8/3; W^4 tt/3; Z 1/2 « 1 (4.10) 

VI. CW EIGENMODES OF A RING BRAGG 
REFLECTOR 


If Bragg-reflecting structure is pumped from one end, 
an incident wave would decay as it propagates by trans¬ 
ferring its energy to a reflected wave due to distributed 
retro-reflection [13]. As one traces both of them in the 
forward direction, one would see the energy of both of 
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them diminishing away from the incidence. However, if 
one pumps the structure from both ends, there could 
be conditions when the combined back-and-forth retro- 
reflection may result in both waves sustaining their en¬ 
ergy without changes. In practical terms, such a double 
pumping could be best attained in a ring (Sagnac) optical 
waveguide resonator [17] coupled to two feeding waveg¬ 
uides delivering counter-propagating pumping. 

Thus another objects of interest here are what we 
call cw ” eigenmodes”, whereby both counterpropagating 
wave do not change as they propagate, i. e. there is no 
energy exchange between them, similarly to the eigen¬ 
modes in and nonlinear wave propagation [18,19]. 
For Bragg reflection, those eigenmodes exist in both lin¬ 
ear and nonlinear case and are allowed only outside the 
Bragg bandgap (which becomes intensity-dependent in 
nonlinear medium). A very interesting aspect of such a 
system is that Bragg eigenmodes should coexist with a 
resonant eigenmodes of a Sagnac resonator, thus creating 
an interesting interplay of resonant effects. 

The eigen-waves propagate with a constant phase ve¬ 
locity, i. e. in eq (1.6) we can assume aj = uje^'^^ where 
the amplitudes Uj, and (unknown yet) normalized per¬ 
turbation of the wave number q are real constants. [Note 
that there is no need to assume time dependence in (5.1), 
since the detuning of the pumping frequency has already 
been taken care for by A in (1.4) or <5 in (1.6).] We will 
be seeking then for a dispersion equation between q and 6 
(which in nonlinear case should also depend on the wave 
intensities uJ) by substituting the solution aj = uje^'^'^ 
into (1.6), and derive that equation as: 

q = V(5-p35'/2)2 _ l-{u\-ul)/2, S = ul + ul{b.l) 
In linear case. S' <C 1, we have 

q = a/(52 - 1 , if (5^ > I (5 2 ) 

The ratio of the amplitudes of both waves is 

uj/uz-j = (-1)V('5 + 3S/2)2 - 1- {6 + 3S/2) (5.3) 

so that in linear case, 

Uj/uz-j = {-ly^/S"^ - 1 - 6. (5.4) 

Both waves here have the same frequency, w = wb(1 + 
fiS), but their wavenumbers are different, kj = ^^[l — 
{—iyqg\. However, in a linear case this doesn’t amount 
to non-reciprocity, since both forward and backward 
waves are coupled, producing a standing-like wave in the 
frame that moves toward the higher intensity wave with 
”k-Doppler” velocity = ujB{ky^— k, qnc/n. 

In particular, at critical points, = 1, we have g = 0, 
hence vop = 0, which indicates a regular standing wave 
with u\ = u\. In a nonlinear case, an intensity induced 
non-reciprocity, the same as in a regular Kerr-line non¬ 
linearity [13] is due to the term u\ — wf in (5.1). 

In a nonlinear case, eq (5.3) suggests a nonlinear con¬ 
nection between the eigen-waves amplitudes, Uj, and the 


frequency detuning of a respective eigen-mode solution. 
Indeed, for any given set of amplitudes uj , the respective 
’’eigen”-detuning is: 

Se^gen = “-S']! -f (miU2)“^]/2 (5.5) 

which, in a linear case, <C 1, is consistent with (5.4), 
written as Seigen = —S/{ 2 uiU 2 ) = —(mi/u 2 + 'U 2 /ui)/ 2 . 
Notice that those eigenmodes in linear case happen only 
outside the Bragg bandgap (5.3); in nonlinear case they 
also determine the boundaries of a new, nonlinearly mod¬ 
ified band whereby and S = 2uiU2' 

5cr = ±l- 3S'/2, (5.6) 

which is red-shifted by 3S'/2 compared to the linear one 
if UK > 0, and blue-shifted if uk < 0. 

VII. DISCUSSION 

The calculations in the previous section are true not 
only for an infinitely long fiber, but also for a ring fiber. 
(A latter case, however, would involve eigen-frequencies 
set up by the ring length.) A ring makes an easy ex¬ 
perimental setup for the observation of those eigenmodes 
by using regular external pumping fibers coupled to the 
ring fiber. Such a ring system and its eigenmodes (both 
linear and nonlinear) may present a considerable interest 
for applications related to Sagnac effect and ring gyro 
based on it, since one my expect a great enhancement of 
the Sagnac effect and hence of the sensitivity of a laser 
fiber-based gyro. 

We didn’t discussed here the stability of those eigen¬ 
modes in strongly nonlinear regime, which can be read¬ 
ily analyzed using small perturbation approach, see e. 
g. [20], similar to the one used in the theory of mod¬ 
ulation instability. It is worth noting however that the 
Bragg nonlinear reflection is known to show optical bista¬ 
bility [21,1,2]; our preliminary research to be published 
elsewhere indicates that above some critical pumping the 
system is actually prone to highly multi-stable and multi- 
hysteretic behavior by forming many quasi-solitons inside 
a finite-length fiber and switching from a N stationary- 
soliton mode to a fV± 1 mode, similar to the propagation 
of strong light in (linearly overdense) plasma layers [10]. 

Furthermore, self-trapped slow or stationary solitons 
in a fiber may be used for the energy storage in com¬ 
puter applications, whereby they can be controlled to 
switch from a stationary state to a high-speed mode and 
used for operational memory or logic operations. For 
that purpose, it may also be of interest to use e. g. 
Erbium, Neodimium, or Ytterium doped amplifying 
fibers or plane waveguides [17] to get them into lasing 
(due to distributed retro-reflection) - and at that slow- 
soliton supporting - mode. 

Another attractive line of further research would be 
to move from ID-fibers to nonlinear 3D-photon crystals, 
which most likely might be instrumental in attaining slow 
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or stationary 3D-solitons via the 3D-Bragg-trapping of 
light and realizing sort of ’’stopped light bullet” that 
would greatly advance and transform the phenomenon of 
a so called ’’light bullet” (we are referring here to a pre¬ 
dicted by Silberberg [22] self-sustained field object mov¬ 
ing with the speed of light of the trapping medium, and 
recently observed experimentally [23]). In this respect, it 
is encouraging that at least a 2D-trapping in the cross- 
section of light beam in the form of random-phase gap 
solitons has been observed in photonic lattices [24], 
Finally, going beyond a regular Bragg-reflection from 
a spatially-periodic structure, one can hypothesize that a 
similar phenomena of slow or stationary solitons might be 
expected in any scattering (but low-absorbing and non¬ 
linear) system, most of all - in those with stochastic scat¬ 
tering. Such a system may have a ” washed-out” bandgap 
edges, yet the main physical factors would remain the 
same - distributed retro-reflection of light forming quasi¬ 
standing wave and facilitating slow-soliton formation via 
nonlinearity. Those properties put it somewhere in be¬ 
tween Bragg-reflector and slightly-overdense plasma [10]. 


Same as above, this idea can be further advanced by 
using amplifying yet scattering media, such as doped ma¬ 
terials. There is a reasonable possibility that a strongly- 
scattering nonlinear system may be able to support 3D- 
self-trapping and to sustain of a long-lived almost unmov¬ 
ing 3D-” hot-ball” reminiscent of a ball-lighting, provided 
there is a sufficient influx of energy from outside to sup¬ 
port inverse population and amplification in the system. 


VIII. CONCLUSION 

A general ’’all-band” intensity profile of slow-moving 
and stationary bandgap (Bragg) solitons for arbitrary 
parameters of of nonlinearity and spacial modulation is 
found to be similar to the profiles of 27r solitons of cascade 
stimulated Raman scattering, including ’’Lorentzian” 
solitons. Nonlinear ”no-energy exchange” eigen-modes 
may become a promising tool in exploration and appli¬ 
cations of such systems. 
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